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ABSTRACT 

We use our Galactic Globular Cluster Catalog (G2C2) photometry for 111 Galactic globu¬ 
lar clusters (GC) in g and z, as well as r and i photometry for a subset of 60 GCs and u 
photometry for 22 GCs, to determine the structural parameters assuming King (1962) models. 

In general, the resulting core radii are in good comparison with the current literature val¬ 
ues. However, our half-light radii are slightly lower than the literature. The concentrations 
(and therefore also the tidal radii) are poorly constrained mostly because of the limited radial 
extent of our imaging. Therefore, we extensively discuss the effects of a limited field-of-view 
on the derived parameters using mosaicked SDSS data, which do not suffer from this restric¬ 
tion. We also illustrate how red giant branch (RGB) stars in cluster cores can stochastically 
induce artificial peaks in the surface brightness profiles. The issues related to these bright stars 
are scrutinised based on both our photometry and simulated clusters. We also examine colour 
gradients and find that the strongest central colour gradients are caused by central RGB stars 
and thus not representative for the cluster light or colour distribution. 

We recover the known relation between the half-light radius and the Galactocentric dis¬ 
tance in the g-hand, but find a lower slope for redder filters. We did not find a correlation 
between the scatter on this relation and other cluster properties. We find tentative evidence for 
a correlation between the half-light radii and the [Fe/H], with metal-poor GCs being larger 
than metal-rich GCs. However, we conclude that this trend is caused by the position of the 
clusters in the Galaxy, with metal-rich clusters being more centrally located. 

Key words: Galaxy; Globular Glusters 


1 INTRODUCTION 

Galactic globular clusters are ancient stellar systems (the major¬ 
ity being formed beyond a ~ 3), containing of ~ 10® stars within 
a volume of ~ 100 pc®. They provide essential information on 
the formation and evolution of the Galaxy (e.g., [Forbes & Bridges| 
| 2010 |l and are natural laboratories for theories of stellar structure 
( jVandenBerg et al.|2013| and references therein). The structure and 
properties of globular clusters bear the imprint of the initial condi¬ 
tions of their formation and interactions with the galactic environ¬ 
ment ( jBrodie & Strader|2006[ ). The high stellar densities (without 

a dark matter halo) make globular clusters invaluable objects for 

* E-mail: Joachimvanderbeke@gmail.com 


the study of N-body dynamics (e.g., [Elson et al.|1987[|Heggie &| 
|Hut|2003[|Trenti et al.|20f0l [Hurley & Shara|2012| l. We also ex¬ 

pect that these complex environments are involved in the formation 
of several stellar exotica, such as blue stragglers ( [Ferraro & Lan-| 
|zoni|200^ jSimunovic & Puzia||2014|(, extreme horizontal branch 
stars ( [Fusi Pecci et al.[1992^, cataclysmic variables and millisecond 
pulsars jBenacquista & Downing | 2013^ , intermediate-mass black 
holes ([Lutzgendorf et al.|2013J and black hole binaries (|Lin et al.| 

|20T^ . 

Structural parameters for globular clusters are needed to ex¬ 
plore correlations between stellar populations, dynamics and the 
Galactic environment. The Milky Way is the only object where 
these questions can be explored in detail, as we can resolve clus¬ 
ters to the level of individual stars on the main sequence and study 
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their internal kinematics jFrank et al.|2012 

Hernandez et aL|2013l 

[Bianchini et al.|20131|Fabricius et al.|2014 

Kacharov et al.|2014|(. 


Most structural parameters for Galactic globular clusters are still 
measured from surface brightness (hereafter SB) profiles derived 
from an inhomogeneous compilation of older CCD and photo¬ 
graphic data I Trager et al.|1995|). The most extensive compilation 


I McLaughlin & van der Marel|2005^ fits these to the classical King| 


1962[)i966| profile for self-gravitatin g systems (|Wilson|1975 and 
other profiles are also used). Recently, [Miocchi et al.| ( |2013^ have 
derived structural parameters based on star counts from Hubble 
Space Telescope (HST) and ground-based photometry for 26 clus¬ 
ters. This latter method is arguably the most accurate (e.g. 

|2003|> but it is observationally expensive (e.g. 


Ferraro 


et al.|1999 


Salinas 


et al.|2012' and cannot be applied to unresolved systems such as 


clusters in all but the nearest galaxies ( |Wang & Ma|2013l >. As long 
as the clusters are bright and well-populated, we expect that SB 
profiles will perform adequately ( [Noyola & Gebhardt|2006| - but, 
see also IGoldsbury et al.||2013) . However, in nearly all instances 
(and almost irrespective of methodology) it is difficult to decide on 
a ’radius’ where globular clusters ’end’. In several cases, clusters 
are seen to contain ’extra-tidal’ stars (e.g., NGC 1851 -|01szewski| 
[eTarillOO^ NGC 5694 - [Correnti et aljllOTT} , while tidal tails 
and other debris are relatively common (e.g., |Grillmair et al.|1995[ 
lOdenkirchen et al.|2001[|Sollima et al.|2011^ . A further complica¬ 
tion is that about 20% of the clusters do not fit smooth King-like 
models, but exhibit a central density enhancement or core-collapse 
( |Cohn|1980) . Given the large ages of these objects, this is expected 
to have occurred in most clusters ([Djorg ovski & Piotto|1993[|Har-| 
|ris|1996^ , but the formation of hard binaries is expected to halt this 
process ^Vesperini & Chemoff]1994[|Fregeau & Rasio|2007^ . 

We have recently completed multi-wavelength observations in 
the Sloan Digital Sky Survey (SDSS) ugriz passbands for the bulk 
of the southern Galactic GCs and supplemented these data with 
northern clusters from the SDSS DR9 ( [York et al.|2000[|Ahn et al.| 
|2012[ l. This results in a total dataset with g and z photometry for 
111 GCs, as well as further r and i photometry for 60 clusters and 
«-band imaging for 22 GCs. We have presented a study of the in¬ 
tegrated griz aperture magnitudes and colours for these objects in 
[ Vanderbeke et al.| (|2014a[ hereafter Paper I) and used these colours 
to improve colour-metallicity relations in | Vanderbeke et aLl ( |2014b| 
hereafter Paper II), the first two papers in this series. Here we use 
the imaging to produce SB profiles and derive structural parameters 
by fitting King models to our homogeneous survey of the Galactic 
globular cluster system. The G2C2 survey has some clear advan¬ 
tages when compared to previous studies: our data are uniform, 
the measurements are based on CCD imaging and the clusters are 
observed during photometric nights ( [Paper I[ l. Moreover, our data 
cover the bulk of the Galactic GCs, were taken in the popular SDSS 
filter system and were analysed carefully and homogeneously. 

The paper is organised as follows. In Section we describe 
the data. Section [^reports on our analysis, presents the stmctural 
parameters and compares those with previous work. At this stage 
we simulate mock clusters to test our fitting algorithm and explore 
possible sources for parameter biases. We discuss the observational 
biases in more detail and use the large images for the clusters from 
the SDSS (which theoretically can be studied to any radius) to study 
the limitations of our dataset by artificially constraining the extent 
of the SDSS SB profiles. We also scrutinise the effects of bright red 
giant stars on the centering of cluster profiles and its consequences 
for the SB profiles and the King model fits, which further confirm 
the findings based on the mock data. In Section |4| we compare our 
[King[([T962l( parameters with studies from the literature using both 


SB and count density profiles. In Section we discuss our results 
and explore their significance for internal colour gradients, stellar 
populations and the interactions of clusters with our Galaxy. Fi¬ 
nally, in Section|^we summarise our conclusions. 


2 DATASET 

The dataset for this paper is based on the Galactic Globular Cluster 
Catalog: a homogeneous imaging survey of a large fraction of the 
Galactic globular cluster system in the SDSS passbands. [Paper I[ 
details the data acquisition process and describes how these data 
were reduced. That study also presented our procedures to correct 
for contamination, how we deal with extinction and calibrate pho¬ 
tometry. We refer to the above paper for details, but give here a 
short summary of the essential information. The bulk of the obser¬ 
vations were carried out on the Cerro Tololo Inter-American Ob¬ 
servatory (CTIO) 0.9 m telescope in a series of observing runs be¬ 
tween 2003 and 2013, with most of the data coming from the 2004 
June season. We performed the standard CCD reduction procedures 
and cosmic ray removal, followed by calibration on the SDSS sys¬ 
tem using stars from [Smith et al.[ ( [2002) . For all clusters we have 
g and z photometry, often taken on several different nights, while 
for an important subset we also have r and i imaging. We supple¬ 
mented this with available clusters within the SDSS footprint (|Ahn[ 
[et al.[[20T2l l. This resulted in a total sample of 111 GCs with gz 
photometry and 60 GCs with ri photometry. 

In [Paper I| we used these data to calculate aperture photome¬ 
try within the (literature) half-light radius (from the 2010 edition 
of the |Harris|[ 1 996[ compilation, which is the version we will use 
throughout this study). This required us to deal with estimating the 
sky level from CTIO images that do not reach to the cluster tidal 
radius. As this is an important issue for the present paper, we dis¬ 
cuss this further below. We also adaptively searched for the cluster 
centre and we similarly describe this in greater detail in this paper. 
Moreover, we experiment in the current study with a new centring 
method based on RGB stars in Section [T3] Finally, we carried out 
removal of contaminating (foreground) stars to obtain clean aper¬ 
ture magnitudes and corrected for extinction using the most recent 
values from [Schlafly & Finkbeiner]j201 1^ . However, in [Paper n| we 
found some evidence that these extinction values may suffer from 
systematic errors, especially near the plane and the Galactic bulge. 

For several GCs in [Paper I| we were unable to carry out aper¬ 
ture photometry because bright red giants in the unresolved clus¬ 
ter centres saturated the CCD. Our purpose in the above publica¬ 
tion was to provide a conservative set of aperture magnitudes and 
colours and we therefore excluded a number of clusters from our 
analysis. Here, we fit King models to the surface brightness profiles 
and we are able to apply an iterative clipping method and profile fit¬ 
ting to derive ’model’ magnitudes. 


3 STRUCTURAL PARAMETERS 
3.1 Surface Brightness Profiles 

We derive radial surface brightness profiles for all clusters and all 
passbands by using conventional aperture photometry, in annuli of 
increasing radius, similar to the procedure of building a curve of 
growth for extended systems. For the CTIO 0.9 m data, where the 
field of view is 13.6' on the side, the annuli are 15 pixels (corre¬ 
sponding to ~ 5.9”) wide, or approximately 4 times the seeing 
disk. We integrated to the edges of the images to obtain the surface 
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brightness profiles to the ~ 6.5' limit allowed by the CTIO 0.9m 
field of view. For SDSS data, we can theoretically integrate the pro¬ 
file to infinity, but we also simulate observations limited to the 4.2' 
and 6.5' field of view used for all other clusters to test the effects 
of the limited aperture on structural parameters. 

We fit the surface brightness profile to a |King| ( [r962| ) model of 
the form: 

^ ^ ([l + (r/r,)2]i/2 " [l + (rt/r,)2]i/2 

where r is the aperture radius, is the core radius and rt the tidal 
radius. The constant k is related to the central surface brightness /i 
as 

1 

[l + (n/r,)2]i/2 

We fit these using a Levenberg-Marquardt non-linear least squares 
algorithm jPress & Schechter|1974[[Bevington & Robinson|1992^ . 
We define a concentration index c' = rt /to simplify our numer¬ 
ical work, while c = log(rt/rc) as in|McLaughlin & van der Marel| 
l |2005^ . We iterated the fit with 5 a clipping about each SB point to 
remove the effects of bright contaminating stars. We also truncate 
the SB profiles, excluding the SB points which are outside the tidal 
radius of the previous iteration. Therefore the outer regions of the 
SB profile were clipped out for some very faint clusters, hence the 
final fit did not rely on the full 6.5' SB profile. The iteration pro¬ 
cess was stopped when the input SB profile vector did not change 
anymore when compared to the one of the previous loop. As a con¬ 
sequence of the 5cr clipping, the central SB point was also ignored 
in some cases, mostly when the SB profile was centred on bright 
RGB stars. 

We show an example King profile for the classical rich cluster 
NGC 104 (47 Tucanae) in Fig.[^ Fits for all other objects are placed 
in an on-line appendix. This rich massive cluster extends much 
further than our limited FOV, thus the concentration index cannot 
be reliably estimated. It is very difficult in any case to constrain 
the tidal radii (and the related concentrations) accurately. Since the 
field of view of CTIO data does not reach to the tidal radius for most 
if not all clusters in the sample, it is inherently difficult to measure 
this quantity. We therefore used an ad hoc estimate, adopting as the 
tidal radius the radius at which the best fit King profiles had flux 
comparable to the mean sky noise. For clusters within SDSS, de¬ 
termining the tidal radius is, in principle, not a problem, since one 
can extrapolate the flux to large distances; even in this case, the 
tidal radius is often difficult to determine jJordan et al.|2009T >. 

We tabulate the derived parameter values based on the CTIO 
data in Table [T] The entire table can be found in the online ap¬ 
pendix. We estimated the errors on these parameters by bootstrap¬ 
ping the derived fits with 100 random points (assuming Poisson er¬ 
rors) and deriving new King model fits to these artificial data. The 
errors in the table are the \a conditional errors on each parameter. 
Where we have more than one observation in one or more filters, 
we always list the King parameters based on the longest one in Ta¬ 
ble or (for similar observing times) the best reduced . Note 
that in some cases (e.g., core-collapsed systems) we were forced to 
exclude the inner annulus from our analysis as it would otherwise 
drive the whole fit because of its small errors: we exclude the cen¬ 
tral SB point when it was more than 1 magnitude brighter than the 
four adjacent SB points. One example is NGC 5927 in the upper 
panel of Fig.|^ This could bias the results for core-collapsed clus¬ 





ters. However, these clusters are by definition not well represented 
by King models in any case. 

We now integrate the King profiles numerically to compute 
a ’model’ magnitude within the half-light radius. These are also 
given in Table [T] For SDSS clusters we give model and aperture 
magnitudes in Table|^ the full table is available online, while an ex¬ 
tract is shown here for guidance. In [P aper i| we presented aperture 
magnitudes based on the half-light from Harris| ( [1996^ . Here we 
determine the half-light radius and we therefore choose to compare 
the model and aperture magnitudes. Generally, these magnitudes 
agree well, with a median magnitude difference (model—aperture) 
of ~ 0.01 mag. 


3.2 Sky values 


In order to follow globular cluster profiles to low surface brightness 
levels, an accurate estimate of the sky value is necessary. Small er¬ 
rors in sky values propagate because of the large areas of the aper¬ 
tures, especially in the outer regions of each cluster. We have used 
the MMM (Mean, Median, Mode) algorithm to measure the sky flux 
in apparently blank regions in the comers of the CCD images: the 
routine is adapted from DAOPHOT ( |Stetson|1987| >, specifically de¬ 
veloped for crowded fields as in our globular clusters. The algo¬ 
rithm clips sky pixels well above the median and combines median, 
mean and mode sky values to obtain a more accurate background 
estimate. In our case, exposure times are relatively short to preserve 
the dynamic range between the bright central regions and the low 
surface brightness wings of the cluster profile; this leads to more 
uncertain sky determinations (see [Paper I| for more details). One 
can compare this with the most commonly used approach to derive 
SB profiles for galaxies jPeng et al.[2002[[2010^ : in some cases the 
outskirts of extended systems are barely above the sky noise and 
even the areal increase in the apertures is offset by the increasing 
noise from the sky, flat field and detector read. This may lead to 
a bias in favour of low effective radii. The problem may be fur¬ 
ther complicated by our assumption of circular symmetry, whereas 
there is evidence of changing ellipticity and position angles in the 
outer regions of some globular clusters ( [Bianchini et al.[2013[ >. Ad¬ 
ditionally, some clusters have stars beyond the tidal radius, while 
others under-fill their tidal region ( [Gieles et al.[ [2010[ |Alexander| 
[& Giele^[2013^ . [Peng et ar| ( [2002[ 20101 argues for the solution 
we have eventually chosen, to estimate the sky independently and 
hold it fixed for the estimate. Finally, in some cases, some very 
bright (sometimes saturated) foreground stars were replaced by the 
sky value. Generally, these stars were already clipped by our fitting 
algorithm. However, the model fits and residuals look nicer when 
masking these saturating foreground stars. It is particularly useful 
to do this for faint sparse clusters, which are sometimes not much 
brighter than the sky itself. 


3.3 Centering errors 

An essential element in our procedure is the choice of an appropri¬ 
ate centroid for the apertures. In [Paper I| we determined the cluster 
centre by using a series of small apertures surrounding the optical 
centre and choosing the position where the flux is maximal (|Noy-[ 
[ola & Gebhardt|2006[ [Bellazzini|2007^ . As a first approximation, 
this is our estimate of the position of the cluster centre. During the 
analysis of the data, it became clear that a few bright stars (mostly 
RGB stars) can strongly affect the derived parameters, mostly be¬ 
cause of their effect on the position of the cluster centre (jGold^ 
[bury et ^[2013^ and especially for poorer systems. Fitting the 
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Figure 1. SB distribution of NGC 104 in the p-band. The left panel focuses on the inner region of the cluster, while the right panel gives details on the outer 
regions. The dashed lines show the 1-sigma deviation of the model, with squares representing the data points that have been clipped iteratively. 
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Figure 2. Surface brightness distribution for NGC 5927 in the z-band. Legend as in Fig.j^ The upper panei shows the fit with naturai (Her) weighting, the 
tower panei exciudes the centrai SB point because it is more than i magnitude brighter than the 4 adjacent SB points. The asterisk for the centrai SB point 
indicates that it was not considered in any of the fitting steps. Boxes indicate the SB points that were not considered in the finai fit. 


centres of GCs accurately has basically been an issue ever since 
King models were fitted to SB profiles. [PJorgovski & Kin^U986^ 
found that about one-fifth of GC cores are not well fit even with 
high-concentration King models. Collapsed cores, which are better 

^ IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astron¬ 
omy under cooperative agreement with the National Science Foundation. 


represented by power-law profiles, are encountered frequently. We 
tested this by creating some mock clusters using IRAhj 

Stars were generated using the|Bahcall & Soneira ^l980|l stel- 
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Table 1. Extract of the King parameter table based on CTIO SB profiles. Central SB uncertainties are pure bootstrapping uncertainties and do not include 
calibration uncertainties, neither the systematic error introduced in . We do not list the King model concentration c, because this parameter could not be 
estimated reliably based on our data. Table [^presents the systematic errors on each of the parameters. Model and aperture (Aper) magnitudes are computed 
within the listed In case of multiple observations, the longest observation was chosen (with exposure time ExpT). The reduced ™d the degrees of 
freedom (DOE) of the fit are also given. SB0i„ indicates if the central SB point was included by the final fit. 


ID 


Mo 

[mag/"2] 

O 

b 

rc 

['] 

cr(rc) 

rh 

['] 


Model 

[mag] 

Aper 

[mag] 

ExpT 

[s] 


DOF 

SBOin 

NGC104 

g 

14.94 

0.00 

0.47 

0.00 

2.62 

0.01 

5.09 

5.08 

60 

0.76 

58 

0 

NGC104 

z 

13.82 

0.00 

0.51 

0.00 

2.52 

0.01 

3.89 

3.85 

60 

1.32 

56 

1 

NGC288 

g 

20.31 

0.00 

1.76 

0.03 

2.19 

0.02 

9.14 

9.12 

270 

0.99 

58 

0 

NGC288 

r 

19.91 

0.01 

2.03 

0.07 

2.05 

0.03 

8.76 

8.70 

60 

0.57 

58 

0 

NGC288 

i 

19.68 

0.01 

2.32 

0.12 

2.01 

0.01 

8.49 

8.41 

60 

0.55 

58 

0 

NGC288 

z 

19.56 

0.01 

2.51 

0.17 

1.88 

0.02 

8.44 

8.39 

60 

0.21 

58 

0 

NGC362 

g 

15.07 

0.00 

0.18 

0.00 

0.97 

0.01 

7.39 

7.35 

60 

0.44 

60 

1 

NGC362 

r 

14.39 

0.00 

0.17 

0.00 

0.80 

0.01 

6.89 

6.92 

60 

0.48 

59 

1 

NGC362 

i 

14.08 

0.00 

0.17 

0.00 

0.83 

0.01 

6.54 

6.58 

60 

0.56 

59 

1 

NGC362 

z 

13.90 

0.00 

0.18 

0.00 

0.70 

0.01 

6.45 

6.50 

60 

0.28 

59 

1 


Table 2. Extract of the GC King parameters and errors for SB profiles based on SDSS data. 


ID 


Mo 

[mag/''2] 

o'(mo) 

rc 

['] 

o-(rc) 

rh 

['] 

(r(rh) 

Model 

[mag] 

Aper 

[mag] 


DOF 

SBOin 

NGC2419 

u 

20.95 

0.01 

0.34 

0.05 

0.60 

0.04 

12.78 

12.76 

0.05 

49 

1 

NGC2419 

g 

19.79 

0.00 

0.36 

0.00 

1.15 

0.00 

10.91 

10.99 

0.76 

52 

1 

NGC2419 

r 

19.33 

0.00 

0.36 

0.00 

1.22 

0.00 

10.41 

10.48 

0.61 

52 

1 

NGC2419 

i 

19.10 

0.00 

0.37 

0.00 

1.19 

0.00 

10.17 

10.23 

0.78 

52 

1 

NGC2419 

z 

19.01 

0.01 

0.37 

0.01 

0.80 

0.01 

10.45 

10.47 

0.08 

50 

1 

NGC5024 

u 

18.89 

0.00 

0.45 

0.00 

0.85 

0.01 

10.11 

10.11 

0.06 

42 

1 

NGC5024 

g 

17.59 

0.00 

0.41 

0.00 

1.07 

0.01 

8.68 

8.68 

0.73 

52 

1 

NGC5024 

r 

17.10 

0.00 

0.38 

0.00 

1.10 

0.01 

8.26 

8.26 

1.23 

53 

1 

NGC5024 

i 

17.19 

0.00 

0.45 

0.00 

1.08 

0.01 

8.16 

8.15 

0.67 

53 

1 

NGC5024 

z 

16.67 

0.00 

0.39 

0.00 

0.92 

0.01 

7.97 

7.96 

0.59 

50 

1 

NGC5053 

u 

23.52 

0.04 

2.32 

0.96 

1.06 

0.01 

13.41 

13.37 

0.01 

56 

1 

NGC5053 

g 

22.43 

0.01 

2.65 

0.31 

2.34 

0.03 

10.91 

10.88 

0.25 

56 

1 

NGC5053 

r 

22.09 

0.01 

2.59 

0.31 

2.24 

0.02 

10.65 

10.65 

0.30 

51 

1 

NGC5053 

i 

21.88 

0.01 

2.79 

0.29 

2.41 

0.05 

10.27 

10.24 

0.36 

56 

1 

NGC5053 

z 

21.81 

0.02 

2.93 

0.87 

1.65 

0.01 

10.80 

10.74 

0.10 

56 

1 


lar luminosity function and distributed over the CTIO FOV accord¬ 
ing to the spatial probability distributions determined by different 
King model density distributions. Fig. shows the distribution of 
the input and c in the IRAF procedure. These input parame¬ 
ters cover the known values from the literature ( |Harris|1996[ l well 
and are extended to somewhat higher concentrations. Because the 
sampling of the clusters can affect the SB profile determination, 
two approaches were used for the stellar densities: images of well- 
populated clusters are two-dimensional projections of clusters with 
10® stars, while more sparse clusters only contain 10^ stars. Ob¬ 
viously, the resulting observed density of the cluster does depend 
on Tc and c itself. No foreground or background contamination is 
included. Exposure times, which determine the signal-to-noise of 
the mock data, were chosen to match our real observations. To de¬ 
rive the cluster centre, we used the above described procedure from 
INoyola & Gebhardt|j2006^ and |Bellazzini|p007) . 

For well-populated clusters measured core radii agree well 
with the input values (Fig. |^. For sparser clusters, the agreement 
is worse, with some showing large offsets. A strong outlier (input 
Tc ~ 1.8', output rc ~ 1') shows how stochastic effects may affect 
our results. We generated this mock GC with 10'* stars, hence it is 
sparsely populated. In fact, the SB profile was centred on a bright 


star, resulting in a strong artificial central peak. This led to a large 
overestimate of the concentration index. If we exclude, arbitrarily, 
the central SB point, this yields more reasonable values. 

In Fig. 1^ we compare the IRAF input concentrations with the 
output of the fitting algorithm. We find that a significant fraction 
of the simulated clusters show large concentration discrepancies, 
especially for sparsely populated clusters. For GCs with tidal radii 
much beyond the FOV, it is impossible to constrain the concentra¬ 
tions, which is not surprising. In general, the output concentrations 
larger than about 2.5 are not reliable. 

Because the distribution of bright stars is stochastic and the 
central surface brightness distributions of King models are flat, er¬ 
rors in the centering may lead to artificially flat profiles ( |Mackey| 
|& Gilmore|2003[ l. Sometimes, however, uncertainties in determin¬ 
ing the cluster center can also produce artificially peaked pro¬ 
files. Therefore, we also experimented with an alternative centring 
method, using the weighted mean position of cluster red giants to 
determine the centre of the cluster. We expect that the red giants 
stars are equally distributed around the centre, hence this method 
should not be biased by the stochastic positions of some bright 
stars. Because this method removes ’artificial’ peaks or dips in the 
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Figure 3. Distribution of the input Vc and c parameters for the IRAF mock 
data. Black asterisks represent our input parameters, small filled grey circles 
show the literature values jHarris|1996). 



Input (orcmin) 


Figure 4. Comparison of the input core radii Vc with the resulting Vc based 
on our fitting algorithm. Well-populated clusters (N = 10^) are indicated 
with a box, other clusters are more sparse (N < 10*^). The dotted line 
shows the one-to-one correspondence. 

SB profile l |Miocchi et al.|2QT3) , it is found that the core radii in¬ 
crease with respect to our original method. 

The following figures compare the core, half-light and central 
surface brightness derived using our original approach and a centre 
based on the positions of red giant stars. The half-light radii com¬ 
pare well between the two approaches, while the central surface 
brightness decreases as artificial peaks in the profile are removed. 

From this, we estimate the systematic error on the parameters 
introduced by the uncertainty in determining the cluster centre and 
list the corresponding values in Table[^ 



Figure 5. Comparison of the IRAF input concentrations with the resulting 
output c based on our fitting algorithm. Legend as in Fig. Due to the 
limited FOV, the concentrations can be overestimated: c > 2.5 are not 
reliable. See text for more details. 



(B07) [orcmin] 

Figure 6 . Comparison of the core radii Vc based on SB profiles obtained 
with the roaming procedure jBellazzini|2Q07| B07) and RGB-based cen¬ 
tres. The red squares indicate the ’red core’ clusters (with ^g—z > 0.1, 
defined in Section [53) . Blue squares indicate the ’blue core’ clusters (with 
Ag _2 < —0.2). The black error bars in the top-left comer show the 2 :-band 
systematic error. See text for more details. 

3.4 Effects of a limited field of view 

We compare clusters in common between the SDSS and our CTIO 
data. In principle, there is no limit to how far a cluster can be fol¬ 
lowed with SDSS photometry and therefore we can test how our 
limited FOV affects model parameters. Note, however, that even 
SDSS data have their limitations: bright RGB stars can saturate the 
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Figure 7. Comparison of the half-light radii rj, based on SB profiles ob¬ 
tained with the roaming procedure jBellazzini|2007| B07) and RGB-based 
centres. Legend as in Fig.[^ See text for more details. 



Figure 8. Comparison of the central SB based on SB profiles obtained with 
the roaming procedure jBellazzini|2Q07| B07) and RGB-based centres. Leg¬ 
end as in Fig.|^ See text for more details. 


Table 3. Systematic errors for the structural parameters due to the stochasti¬ 
cally distributed RGB stars and the uncertainty on the centre determination. 
For the ri parameters, the systematic uncertainties for the 2 :-band should be 
adopted for all parameters. See text for more details. 


Filter 

'^syst (A^o) 

^systiXc^ 

^ systiXh) 


[mag/"2] 

['] 

['] 

g 

0.228 

0.077 

0.101 

z 

0.289 

0.118 

0.191 


SDSS CCD; moreover, the mosaicking of the SDSS data can result 
in artificial sky gradients. 

We use the SDSS mosaics to determine the surface bright¬ 
ness profile and fit this to a King model. We then artificially restrict 
the field to radii of 4.2' and 6.5' (the total CTIO FOV), to esti¬ 
mate the influence of the field of view on the model fits. In general, 
we find that the smaller field of view biases the results towards a 
smaller core radius, but using the full 6.5' FOV solves this issue 
(Fig. 1^. Note that tidal tails were found in NGC 5053 (|Jordi &| 
|Grebdl2010| l. 

In Fig.[T^we present a comparison of the half light radii based 
on the full mosaicked SDSS frame and the artificially truncated 
profiles. These two quantities compare well, which suggests that 
the half-light radius is also generally well determined. The scatter 
on the one-to-one correspondence is again reduced when using the 
6.5' SB profile instead of the 4.2' SB profile. 

Fig. [m shows a comparison of the structural parameters for 
clusters with both SDSS (full images) and CTIO data. The central 
surface brightnesses are in good agreement (with one exception - 
NGC 7078 - whose profile is notoriously difficult to fit -|Newell &| 
|Oneil|1978p . Core radii are also in good agreement, while concen¬ 
trations from our data are generally found to be unreliable. 


3.5 Comparison of observations in the gz filters 

Different filters are sensitive to different stars. However, if the King 
parameters truly reflect the star count densities, we do not expect 
substantial differences between the structural parameters based on 
observations in the blue or red. 

Fig. [T^ compares the core radii and half-light radii based on 
the g and 2 filters. We do not find evidence for a systematic trend 
in the core radius or the half-light radius. 

GCs close to the Galactic plane are known to suffer from dif¬ 
ferential reddening ( [Alonso-Garcfa et al.|201 1| >. Observations taken 
at longer wavelengths are less affected by the extinction. Therefore, 
we expect a correlation between the differential reddening and the 
gz structural parameter differences, if the differential reddening 
plays a significant role in the determination of the King parame¬ 
ters. However, we do not find any correlation and conclude that the 
differential reddening does not strongly bias our parameter deter¬ 
minations. 


4 COMPARISON WITH THE LITERATURE 


We now compare our results with previous determinations of struc¬ 
tural parameters from the literature. First, we compare with the 
20 10 edition of | Harris (|1996i, which is largely a compilation of 
the |McLaughlin & van der Marelj(|2005|l SB profi le study. Then we 
also compare with the results of |Miocchi et ^ ( |2013| >, who used 
star count density profiles to determine the structural parameters. 

It is clear that some clusters have lower quality data, resulting 
in more uncertain parameters. We therefore restrict our analysis to 
a smaller, ’’clean” sample where we can trust our determinations of 
the main structural parameters. We impose following constraints: 


• Tc > o'(rc): The core radius has to be significantly different 
from zero (based on the bootstrap errors, not including the system¬ 
atic errors). 

• go <20 mag/"^: clusters fainter than this SB have very large 
photometric errors. SB points are often not significantly different 
from zero. 




































8 J. Vanderbeke et al. 


Limited FOV 4.2 arcmin Limited FOV 6.5 eremin 




Figure 9. Comparison of the core radii Vc based on the full mosaicked FOV and artificially limited SB profiles of 4.2' and 6.5' (left and right panel, 
respectively) for all SDSS GCs belonging to the clean sample introduced in Section|^(e.g. 19 GCs for g, 18 GCs for z). The artificial cuts limiting the SB 
profile radius are indicated with the dashed lines. The dotted line represents the one-to-one correspondence. It is clear that the scatter is significantly reduced 
using the 6.5' SB profiles, stressing the importance of using the entire CTIO FOV. See text for more details. 

Limited FOV 4.2 arcmin Limited FOV 6.5 arcmin 




limited FOV (arcmin) limited FOV (arcmin) 

Figure 10. Comparison of the half light radii rf^ based on the full mosaicked FOV and an artificial limited SB profile for the SDSS GCs belonging to the 
clean sample introduced in Section|^(e.g. 19 GCs for g, 18 GCs for z). The dotted line represents the one-to-one correspondence. It is clear that the scatter is 
reduced using the 6.5' SB profiles, stressing the importance of using the entire CTIO FOV. See text for more details. 


These conditions are met for 12, (80, 41, 42, 87) in the u 
(g, r, i, z, respectively)-band. We have experimented with SB cri¬ 
teria including the effects of foreground extinction (basically im¬ 
posing magnitude limits on the apparent SBs, not corrected for ex¬ 
tinction), but this did not strongly affect the selected clusters for 
the clean sample. Therefore, we do not introduce restrictions on 
the foreground extinction as this would automatically exclude all 


GCs close to the Galactic plane. However, some disk/bulge clusters 
were removed from the sample, either because of the strong stellar 
contamination or because of the low apparent SB points (often not 
significantly higher than the sky, partly due to the foreground ex¬ 
tinction). Because the number of clusters is largest in the grz-bands, 
we will focus on these filters during the analysis and discussion. 
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Figure 11. Comparison of the King model parameters (central SBs fig, core radii rc and concentrations c) based on the full mosaicked SDSS FOV and on the 
CTIO data. The dotted line represents a one-to-one correspondence. The z-band systematic errors are shown by black error bars. 
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Figure 12. Comparison of the King parameters derived based on observations with the g and z filters. The black error bars present the corresponding systematic 
errors, as listed in Table]^ 


4.1 Comparison to parameters based on SB profiles 

Fig. [^compares c and to the literature ( |Harris|1996^ . The 
left panel of the figure compares the core radii (based on the g- and 
2 -bands). The match for the well-populated clusters is very good, 
with the exception of NGC 5139. For this cluster, our fitting algo¬ 
rithm clipped the central SB point (for both g and 2 filters), which 
accounts for the larger value of the core radius. However, for poor 
clusters (indicated in the figure), our resulting can be biased. The 
middle panel compares our concentrations with |Harns] ( | 1996^ . For 
the bulk of our sample, the concentrations are unreliable and will 
therefore not be discussed further. 

In the right panel of Fig.[T^we show our half-light radii vs. the 
values in the |Harris|jl996| compilation. For smaller clusters (rj, < 
1'), the data fall close to the 45-degree line, albeit with somewhat 
large scatter. For clusters subtending a larger angle on the sky, we 
find that our rn are systematically lower, probably because of our 
poorer value of rt or a poor sky determination (see Section [0] >. 

As an additional check, we selected a representative subsam¬ 
ple of CTIO clusters and fitted Sersic profiles in pixel space with 
BUDDA ( |de Souza et al.|2004}|Gadotti|2008^ . Again, we found that 
the effective Sersic radius (equivalent to the radius encompassing 
half of the cluster light) was significantly smaller (about a factor 
three) than the literature values. Both approaches suffer from the 
sky determination issues related to the limited FOV of the CTIO 
data. Nevertheless, the centre determination is incorporated in the 


BUDDA fitting algorithm, hence we do not expect a strong bias due 
to central RGB stars. 

The King model parameters are fitted simultaneously. There¬ 
fore the discrepancies between the resulting parameters and the val¬ 
ues given by |Harris] ( |1996^ could be correlated. In Fig.[T^we study 
these correlations in more detail and conclude that no correlations 
between Arc, Arji and Ac are present. 

4.2 Comparison to star count density profiles 

Star count density profiles mainly trace the abundant main- 
sequence stars in GCs. SB profiles are vastly affected by bright 
evolved RGB stars, that sink to the cluster centre due to mass seg¬ 
regation jFregeau et al.|20d2} . In contrast to SB profiles, star count 
profiles are not biased by the presence of sparse, bright stars ( |Noy-| 
|ola & Gebhardt|2006| l. Therefore, the latter are believed to be the 
most reliable way to determine structural parameters. However, our 
CTIO and SDSS data does not have the required resolution to de¬ 
rive such radial stellar density profiles. Based on a combination of 
high-resolution HST observations and ground-based observations, 
|Miocchi et al.| ( [20T3) ) were able to construct star count profiles and 
derived structural parameters for 26 Galactic GCs. 

Fig. presents a comparison between our King model pa¬ 
rameters and the parameters derived by |Miocchi et al.| ( |2013^ . Our 
concentrations, shown in the left panels, again prove unreliable. 
The middle panels compare our newly derived core radii with the 
















































10 


J. Vanderbeke et al. 


NGC51^§ 


■ NGC4p, 
NGC5466 

9 




0 12 3 4 

(arcmin) Harris 1996 



c Harris 1996 



0 1 2 3 4 5 

(arcmin) Harris 1996 


Figure 13. Comparison of the structural parameters between our results (for the g- and ^-band) and the literature jHarris|1996) . Filled circles are values based 
on CTIO SB profiles, supplemented with SDSS values represented by asterisks. The dashed line indicates the one-to-one correspondence. High-concentration 
clusters (with c > 2.5) are indicated with red boxes in the left and right panels. Concentrations higher than ~ 2.5 are unreliable, which was also concluded 
from the simulations in Fig.[^ Therefore, clusters with c > 3 are not included in the concentration panel. However, both and rj, are recovered well for the 
bulk of the cluster sample. The systematic error (given in Table[^ is illustrated by the black error bar in the top-right corner of each panel. See text for more 
details. 


core radii derived by |Miocchi et al.| ( f2013^ based on King model fits 
to star count density profiles. In general, these compare well, with 
some exceptions which are indicated in the figure. The right panels 
of Fig. [T^ show a comparison between our King model rn and the 
effective radii derived by |Miocchi et al.| j2013| >. The scatter is 
rather large with some significant outliers. 


5 DISCUSSION 

5.1 Correlations between size and metallicity 


Several extragalactic studies found that red (metal-rich) GCs are 
systematically smaller (about 20%) than blue (metal-poor) GCs 
(see, e.g., |Kundu & Whitmore|1998[|Kundu et al.|1999[|Puzia et al.| 


|1999[ [Larsen et al.|200lj r Larsen & Brodie ( 2003| > argue that the 

difference is originated by a projection effect, due to a correlation 
between GC size and Galactocentric distance, combined with dif¬ 
ferent radial distributions of the metal-poor and metal-rich GC sub¬ 
populations (with the metal-rich GCs being more centrally located; 


see e.g. Kissler-Patig et al.|1997[|Lee et al.j 1998][Cote et al.|200l| 


[Dirsch et al.|2003 l. At larger distances from the galactic centre, the 
projection effect is less strong, hence the trend should disappear 
in the galactic outskirts. However, |Harris| ( f2009| l studied GCs from 
massive galaxies and found that the size differences between red 
and blue GCs remained at large distances. An alternative scenario 
was proposed by |Jordan| ( |2004| l, who advocate that the correlation 
between colour (metallicity) and is a consequence of mass seg¬ 


regation and the longer lifetimes (for a given mass) of more metal- 
poor stars. 

In Fig. |I6| we show the half-light radii as a function of the 
metallicity |Hams TT996t, colour-coded by age (which are largely 
obtained from Forbes & Bridges|2010| with some other references 
as described in Paper II[ . The correlation between both param¬ 
eters is rather weak (Spearman’s rank correlation coefficient is 
Ps,g ~ —0.4 with a significance level of 4.24 x for the g- 
band and ~ —0.23 with a significance level of 3.12 x 10~^ 
for the a-band). Nevertheless, [Jordan] ^2004] ! stress that large sam¬ 
ples are needed to accurately determine the average behaviour. Al¬ 
though Galactic GCs are not so numerous, we still try to recover 
the average behaviour by showing the median [Fe/H] and ru in 
0.5 dex metallicity bins. These median values largely agree with 
the trend predicted by [Jordan] ( [2004^ : rh,g ~4.1 pc, 3.1 pc, 2.5 
pc, 1.9 pc and 2.1 pc for the different metallicity bins {rh,z ~3.8 
pc, 2.9 pc, 2.8 pc, 2.6 pc and 2.7 pc), thus the difference be¬ 
tween g and 2 -band for the different metallicity bins is 0.3 pc, 0.2 
pc, —0.3 pc, —0.7 pc and —0.6 pc. However, the actual size dif¬ 
ferences are much larger than predicted by [Jordan| ( [2004[ . The in¬ 
terpretation of the trend between cluster size and iron abundance 
is not that straightforward. Therefore, we will revisit this issue in 
Section [53] where we discuss the correlation between the GC size 
and the Galactocentric distance. 

Note that for the p-band, the youngest (i.e. < 8 Gyr) clusters 
are systematically below the median values for the corresponding 
metallicity bin, while this is less clear for the 2 -band. 
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Figure 14. Coirelations between the differences of the King parameter determinations when compared to the literature (Ap = pthis work — PHanis 1996 for 

p = rc,c, rh). 


5.2 Correlations between age and core radius 

[Mackey et al.| ( |2008| l found a trend between age and core radii 
for the GCs in the Small and Large Magellanic Clouds, with the 
young clusters (with ages between 10 Myr and 1 Gyr) showing a 
much smaller spread than the old GCs (older than 1 Gyr). These 
authors investigated two physical processes leading to large-scale 
core expansion: mass-loss due to rapid stellar evolution in a primor- 
dially mass-segregated cluster and dynamical heating due to stellar 
mass black holes. 

Galactic GCs are preferentially old jForbes & Bridges|2010|l, 
hence one would expect to see a large range of rg. However, Fig.| l7| 
shows that the majority of the clusters have small g-band core radii. 
[Mackey et al.|(2008| > believe this is a consequence of the positions 
of the GCs in the Galaxy: because a large fraction of the MW GCs 
reside in the inner Galaxy, tidal forces are expected to rapidly de¬ 
stroy loosely bound clusters. The same conclusions can be drawn 
based on the 2 -band core radii. Both NGC 2419 and NGC 5139 
(Omega Cen) have properties that hint to an extragalactic origin, 
which might explain their offset to the general trend. 


5.3 Correlations with the Galactocentric distance 

[van den Bergh et'SII ( [1991^ and [McLaughlin[ ( [2000^ have pro¬ 
posed the existence of a relation between the Galactocentric dis¬ 
tance (Roc) and the effective radius, arising from the tidal trun¬ 
cation (but see also |Kundu et al.[ ( [1999[ >). Recently [Ernst & Just[ 
( [2013^ provide a more theoretical approach on the implications of 
the correlation between the half-mass and the galactocentric radius. 
[Miocchi et al.|([20T^ recover this correlation and show that it does 
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not depend on other cluster properties, confirming its likely dynam¬ 
ical origin. [Puzia et al.[ ( |2014^ studied the GC system of NGC 1399 
and pointed at the influence of the GC orbit distribution on the evo¬ 
lution of the structural parameters. Fig.[^shows the observed cor¬ 
relation between the half-light radius vs. Galactocentric distance in 
our g- and 2 -band data. Based on a robust fit to the data points, we 
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Figure 15. Comparison between our King model parameters (c, Vc, Vh) and the parameters based on King model fits to star count density profiles ^Miocchi] 
|et al.|2Q13) . The dashed line indicates the one-to-one correspondence. Concentrations higher than ~ 2.5 are unreliable, which was also concluded based on 
the mock data in Fig.j^and from Fig.|13|The black enor bar in the top-left corner of each panel illustrates the systematic error given in Table|^ 


find: 

logrh = -1.30 ± 0.10 + (0.45 ± 0.05) x logi?Gc (2) 

for the (/-band (Spearmans rank correlation coefficient ps,g is about 
0.68 with a significance level of 2.44 x 10“^^), which is in good 
agreement with the empirical power law found by |van den Bergh| 
|et al.|(T99T} and|Mackey & van den Bergh|(|2005ll, and, in reason¬ 
able agreement with the scaling relation found by [Miocchi et al.| 
( |2013^ . The latter study, though limited in sample size, covers sev¬ 
eral GCs at large Roc- 

Surprisingly, for the z-band, we find: 

log Vh = -0.65 ± 0.09 + (0.29 ± 0.04) x log Rgc, (3) 

which is significantly different (ps,z ~ 0.53, significance level 
~ 1.49 X 10“^). To further investigate the origin of the offset, 
we made similar robust fits to the clean r and i-band samples. We 
obtain: 

logr,, = -0.98 ± 0.13 + (0.38 ± 0.07) x logi?Gc (4) 

for the r-band (ps,r ~ 0.57, significance level 7.71 x 10“®) and 

logr,, = -0.76 ± 0.12 + (0.33 ± 0.06) x logi?Gc (5) 

for the i-band (ps,i ~ 0.60, significance level 2.69 x 10“®), hence 
we find evidence for a decreasing slope for redder bandpasses. 
Clusters are redder outwards, either because of tidal truncation or 
because shocks with the Galactic centre environment speed up mass 
segregation. Colour gradients will be discussed in full detail in Sec¬ 
tion [5^ It will become clear that the colour gradients are mainly 


caused by the influence of RGB stars and are not representative for 
the cluster count densities. 

[Georgiev & B6ker| ( |2014| l find that nuclear star clusters (with 
sizes ranging between GC and UCD sizes) have smaller effective 
radii when measured in bluer filters. Now we will study the size 
differences of GCs when measured in different filters, and, we will 
try to link those differences to the discrepancy in the Rgc — 
relations. 

In Section l5Tl we found tentative evidence for the existence 
of a median size difference related to the cluster metallicity: metal- 
rich clusters, which are preferentially located close to the Galactic 
centre, have larger Vh when measured in redder filters. Here we 
define Arh,g-z as the difference (in arc minutes) between the g- 
band and 2 -band half-light radii. For the entire sample, the median 
(mean) of Arh,g-z is about —0.06' (—0.07', respectively), which 
is consistent with a zero difference within the errors. However, for 
GCs with Rgc lower than 5 kpc, the median (mean) of Ar^^g-z is 
about —0.19' (—0.23', respectively), while for GCs at larger dis¬ 
tances to the Galactic centre, the median (mean) of Ar^^g-z is 
about 0.09' (0.08', respectively). Therefore, close to the Galactic 
centre, clusters have redder outskirts (or, equivalently, bluer cen¬ 
tres). For GCs with E{B — V)< 0.1 {E(B — V)> 0.1), the me¬ 
dian of Arh,g-z is about 0.08' (—0.15', respectively), hence GCs 
at low reddening appear larger in g than in z, while GCs at high 
reddening appear smaller in g than in 2 . Similar correlations are 
found when using the absolute distance above the Galactic plane 
instead of the extinction value. 

To further scrutinise the effects of metallicity and tidal disrup- 
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Figure 16. Correlation between and [Fe/H]. The black filled stars present the median values for metallicity bins of 0.5 dex. 
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Figure 18. Correlation between the half-light radius r^, and the Galactocentric distance Rqc for the g- and 2 :-band. Small green filled diamonds indicate the 
GCs with [Fe/H] > —1. The solid lines present robust fits to the data points of the clean sample (given by Eqs. [^and|^. 


tion on the cluster size, we make separate fits for the metal-rich and 
metal-poor suhsamples of the clean sample. Because the metal-rich 
GCs do not span a wide range in Galactrocentric distances, we fix 
the slope of the relation to the value found in Eqs. 1^ and and find 
that the intercept for metal-rich and metal-poor subsamples is not 
significantly different. This points to a scenario in which the origin 
of the size difference is related to the Galactocentric distance, and, 
not to [Fe/H] (hence in contrast to the tentative evidence found in 
Fig. p). 

Puzia et ar] (|2014|) found a flattening of the Rgc — ru corre¬ 
lation in NGC 1399 beyond ~ 20kpc, while [Miocchi et aL|j2013^ 
did not find evidence for a flattening of the latter relation in the out¬ 
skirts of the Milky Way. Some GCs at ~ 20 kpc with good fits are 
located below the general trend and could point at a flattening of 
the correlation. Nevertheless, the scatter on the correlation is sig¬ 
nificant. 


To detect correlations between the deviation to the general 
trend and other parameters, we define 

Dg = -1.30 -I- 0.45 X log(iZGc) - log(rh) (6) 

for the g-band (and a similar definition for the a-band; see |van d^ 
|Bergh|2012| for a similar approach). We searched for correlations 
between this new parameter and age, HB index, absolute magnitude 
and [Fe/H], but did not find any significant trends and confirm 
therefore the purely dynamical origin of the correlation between 
rh and Rgc, concurring with the recent literature jvan den Bergh| 
|2012[|Miocchi et al.|2013| >. We also looked for correlations between 
the absolute distance D and the absolute distance to the Galactic 
plane, which could reveal some insights on the influence of tidal 
shocks during disk passages on the GC size. However, we did not 
find any clear trends. 

Fig. [T^ shows the correlation between the absolute distance 
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Figure 20. The 3 -band (grey solid line) and z-band (black dashed line) 
distributions of the core radii. 



Figure 21. The gz distributions of the half-light radii. Legend as in Fig.|20| 


above the Galactic plane and the half-light radii. The clusters lo¬ 
cated closer to the disk are generally smaller than those well above 
the Galactic plane, which can be ascribed to tidal stripping. For the 
3 -band, the Spearmans rank correlation coefficient ps,g is about 
0.61 with a signihcance level of 2.66 x 10“®. For the 2 -band, the 
correlation is less strong: for the entire sample we hnd ps,z ~ 0.46 
with a signihcance level of 8.66 x 10 “®. 


5.4 Distributions of the structural parameters 


The structural parameters of a cluster can be taken as a broad mea¬ 
sure of its dynamical evolution. In |Trager et al.| ( [T995^ the distribu¬ 
tion of core radii of globular clusters was found to be bimodal, with 
about 20 % of clusters showing central light excesses consistent 
with core collapse. Cusps and similar excess light features appear 


to be relatively common in globular clusters in our Galaxy 1 Ves- 


perini & Trenti|| 2010^ and the Large Magellanic Cloud ( [Mackey 
& Gilmore|2003 1 . In Fig. we show the gz distributions. We 


do not hnd clear signs for bimodality, in contrast to|Trager et al.| 

( fT9^ . 

Fig. [^presents the gz distributions of the half-light radii. We 
hnd a median of 2.7 pc (2.9 pc) for the 3 -band ( 2 -band, respec¬ 
tively), which is in good agreement with rji ~ 3 pc found for extra- 
galactic GCs in the Virgo and Fornax clusters of galaxiesJJoidm 


et al. 


20141. 


2005[|Masters et al.|2010[ see also e.g. Fig. 16 ofjPuzia et al. 


Standard dynamical models of GCs make clear predictions 
of the evolution of the ratio between the core and the half-light 
radii (see [Trenti et al.||2010l [Miocchi et al.||20T^ and references 
therein). Fig. [^presents the distribution of the ratio between the 
core and the half-light radii. We hnd a peaked distribution centred 
at Tc/vh ~ 0.4, but do not recover the bimodal distribution found 
by |Miocchi et al.H2013[ l, with a peak in r^jrh at about 0.3. Both 
our and their rdrh values are in agreement with expectations from 
simulations of cluster dynamical evolution. 



Figure 22. Distribution of the ratio between the core and the half-light radii. 
Legend as in Fig.|20| 


5.5 Colour Gradients 

In galaxies, colour gradients are interpreted as a metallicity gradi¬ 
ent jTamura & Ohta|20M||La Barbera et al.|2010^ . If colour gradi¬ 
ents exist in GCs, these would not be linked with the metallicity, as 
these objects have a largely homogenous metallicity (although vari¬ 
ations in light element abundances are omnipresent). However, if 
dynamical processes affect stellar populations, for instance through 
the creation of blue stragglers in cluster cores ( [Ferraro & Lan-j 
[zoni[200^ or strippin g of red giants to produce AGB-manque stars 
( [Pasquato et al.[2013[ l, the formation of cataclysmic variables, mil¬ 
lisecond pulsars and intermediate mass black holes, we might be 
able to detect the presence of these objects through colour vari¬ 
ations as a function of radius. [Djorgovski & Piotto[ ([1993[( found 
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Figure 19. Correlation between the absolute distance above the Galactic plane and . 


that post core-collapse clusters are bluer in their centre and there 
were no colour gradients in the opposite direction. |Djorgovski et al.| 
([T 99 TJ found no blue cores in clusters with flat inner prohles, con¬ 
firming the correlation between internal dynamics and stellar pop¬ 
ulations. This was attributed to the stripping of red giants in cluster 
cores to produce extended HBs ( [Pasquato et al.|2013| l and the for¬ 
mation of blue stragglers by stellar collisions. 

However, |Sohn et ^ ( |1996[ ) found colour gradients even in 
clusters that fitted the King profile and suggested that internal dy¬ 
namics may produce extended HBs and lead to colour gradients. 
|Sohn et ^ jl998| l also found both red and blue cores in both nor¬ 
mal and post core-collapse clusters. For example, unlike [^or gov-| 
|ski et al.| ( |1991t , |Sohn et al)] ( |1998| l detect a red core in NGC 2808. 

We find general flatness of the gradients beyond the half-light 
radius, which is somehow expected as we do not expect strong stel¬ 
lar population gradients. However, stronger gradients are apparent 
closer to the cluster centre. To study the colour of the cluster cores 
and search for correlations with other GC parameters, we define a 
new parameter: 

Ag —2 = {g ~ z)r^ ~ {9 ~ (7) 

with {g — z)rc the g — z colour at the average gz core radius (similar 
for the half-light radius). This parameter is negative for blue cores 
and positive for red cores. 

In Fig.|^we present the distribution of this newly defined pa¬ 
rameter, which is rather symmetric. Clusters with extremely pos¬ 
itive values (Aj ,_2 > 0.3) are NGC4833 and NGC6584, ex¬ 
tremely negative values (Ag-z < —0.4) are found for NGC6235, 
NGC6352 and NGC6681. These clusters are close to the Galactic 
disk or bulge, which results in SB profiles polluted by Milky Way 
stars and complicates the determination of the structural parame¬ 
ters. 

Clusters in the "blue core tail” (with —0.4 < Ag-z < —0.2) 
are NGC 6121, NGC 6342, NGC 6453, NGC 6535, NGC 6638, 
NGC 6637, NGC 6642, NGC 6712, NGC 6864 and NGC 6981. 
Clusters in the ’’red core tail” (with 0.1 < Ag-z < 0.3) 
are NGC 2298, NGC 5634, NGC 6171, NGC 6273, NGC 6341, 
NGC 6760 and NGC 6779. For a fraction of these clusters, we also 
determined the King parameters based on a SB profile centred with 
the RGB method, a procedure which should not be susceptible to¬ 
wards centrally located bright stars (see Section [T^. About one 



Figure 23. Distribution of Ag—z (defined in Eq.|7j. The dashed line repre¬ 
sents the median value of Ag-z- See text for more details. 


third of the ’’tail” clusters are located towards the Galactic bulge, 
prohibiting an RGB-based centre determination. For a couple of 
clusters, only poor colour-magnitude diagrams (CMDs) were ob¬ 
tained and a centre determination based on the RGB was not possi¬ 
ble. 

It is generally difficult to study colour gradients, especially 
close to the cluster centre, because only a couple of stochastically 
distributed RGB stars can alter the SB profiles and originate colour 
gradients. This is especially true for low-density clusters, but can 
also affect clusters as bright as NGC 104: RGB stars and poor cen¬ 
tring can affect the resulting King parameters drastically. 


6 SUMMARY 

As a part of the G2C2 project, we use in the current study our op¬ 
tical dataset (presented in [Paper I]) to derive structural parameters 
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fitting [Km^ ( | 1 962| ) models to the GC SB profiles. We present struc¬ 
tural parameters in the g and « filters for 111 Galactic GCs, while 
we also include r and i-band values for 60 clusters and M-band es¬ 
timates for 22 GCs. For some clusters it was not possible to fit a 
representative King model to the SB profiles. Neglecting the un¬ 
realistic fits (judged by visual inspection), we present u, g, r, i 
and 2 structural parameters for 18, 94, 50, 48 and 95 GCs, respec¬ 
tively. Nevertheless, in order to maintain the focus on the GCs with 
soundly determined parameters in the discussion, we find it nec¬ 
essary to introduce some parameter criteria in the discussion (Sec¬ 
tion]^. For the u (g, r, i, 2 ;)-band, the stmctural parameters of 12, 
(80, 41, 42, 87) GCs satisfy these standards. 

Because the bulk of our dataset was obtained with CTIO 0.9 m 
telescope observations, suffering from a limited FOV, we exten¬ 
sively discuss the effects of a limited field-of-view on the derived 
King model parameters. In general, the resulting core radii are in 
good comparison with the current literature values. However, our 
half-light radii are slightly underestimated when compared to the 
literature. The concentrations (and therefore also the tidal radii) are 
poorly constrained, partly due to the limited radial extent of our SB 
profiles. Moreover, RGB stars are biasing the centring procedure 
and can provoke strong central SB cusps, which further contribute 
to the overestimation of the concentration parameter. The issues 
related to these bright stars are scrutinised based on our photomet¬ 
ric data and simulated clusters. Logically, the effects of the ran¬ 
domly distributed RGB stars on the fitted King models are stronger 
in sparsely populated clusters, for which a star count density ap¬ 
proach is recommended. Colour gradients are examined and can 
also be related to RGB stars. 

We recover the known relation between the half-light radius 
and the Galactocentric distance for the g-band, but find a lower 
slope for redder filters. We did not find a correlation between the 
scatter on this relation and other cluster properties. We find ten¬ 
tative evidence for a correlation between the half-light radii and 
[Fe/H], with metal-poor GCs being larger than metal-rich GCs. 
However, it turns out that this trend is due to the relation between 
the half-light radius and the Galactocentric distance, with metal- 
rich clusters being more centrally located than metal-poor clusters. 
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